function [ out ] = Igain( pred, values, problemID )
%UNTITLED6 Summary of this function goes here
%   Detailed explanation goes here
N = problemID;
f = @(x) pred(x);
id = f(values);
n = sum(id(logical(N))==1);
p = sum(id) - n;
M = max(size(values));
K = sum(N);
out = 1/hygepdf(n, M, K, sum(id));
end

